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Abstract 

External loads typically have a indirect influence on phonon curves, i.e., they influence the 
phonon curves by changing the state about which linearization is performed. In this paper, we 
show that in nanotubes, the axial load has a direct first-order influence on the long- wavelength be- 
havior of the transverse acoustic (TA) mode. In particular, when the tube is force-free the TA mode 
frequencies vary quadratically with wave number and have curvature (second derivative) propor- 
tional to the square-root of the nanotube's bending stiffness. When the tube has non-zero external 
force, the TA mode frequencies vary linearly with wave number and have slope proportional to the 
square-root of the axial force. Therefore, the TA phonon curves - and associated transport prop- 
erties - are not material properties but rather can be directly tuned by external loads. In addition, 
we show that the out-of-plane shear deformation does not contribute to this mode and the unusual 
properties of the TA mode are exclusively due to bending. Our calculations consist of 3 parts: first, 
we use a linear chain of atoms as an illustrative example that can be solved in close-form; sec- 
ond, we use our recently-developed symmetry-adapted phonon analysis method to present direct 
numerical evidence; and finally, we present a simple mechanical model that captures the essential 
physics of the geometric nonlinearity in slender nanotubes that couples the axial load directly to 
the phonon curves. We also compute the Density of States and show the significant effect of the 
external load. 



1 Introduction 

The low-dimensionality of nanotubes causes unusual features in the phonon curves due to geometric 
effects. In this paper, we examine particularly the transverse acoustic (TA) phonon mode and its slope 
and curvature (second derivative) in the long wavelength limit. 

Previous works have reported four acoustic modes of vibration for SWNT. There is consensus that 
two of the modes (longitudinal and twisting modes) have linear dispersion in the long wavelength 
limit. However, there has been discussion about the form of the other two TA modes which are doubly 
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degenerate. E.g., several papers [!Il|2l|3l|4l[5l[6l|71[8l|9l report linear dispersion for the TA modes. 
These authors do not analytically show that the TA modes should be linear, but argue that nanotubes do 
not have modes analogous to the out-of-plane "bending" mode of graphene, and thus, nanotubes must 
have linear long wavelength acoustic dispersion behavior. Among these works, two recent papers EIH 
state that quadratic dependence is found in some other papers, but argue that this may be attributed to 
numerical errors. 

In contrast, more recent papers flOj [TTl [T2l [T3l [T4l [T5l [T6l [TTll report a quadratic dependence 
of frequency on the wave number. The results presented in these papers are obtained by applying 
lattice dynamics methods to either tight-binding, ab-initio or empirical potential models. In still other 
papers [[TSl [T9l l20l |2T]| analytical methods are used to argue that the TA mode should be quadratic. 
In particular, the papers ifTSl |T9l l20ll make an analogy between the phonon dispersion relation and the 
bending energy or sound velocity in a continuous thin tube. The paper [21 J develops a spring model 
and performs lengthy calculations for zig-zag and armchair nanotube^ Among these papers, the 
works [fm [T3l [T6l l20l IT9l [T8l also mention that others have observed linear dependence of the TA 
dispersion relations. However, these papers do not provide a reason for the discrepancy but rather only 
postulate that the use of unsuitable force constant parameters may lead to linear TA phonon dispersion. 

In this paper, we examine the question of the TA phonon dispersion through a combination of 
explicit calculations in model systems as well as numerical calculations. Our model system consists 
of a linear chain of atoms as an illustrative example that can be solved in close-form. The numerical 
calculations are based on recently-developed symmetry-adapted methods for nanotubes that allow the 
examination of bending modes [|22l l23l . Our calculations show that in nanotubes, the axial load has 
a direct first-order influence on the long-wavelength behavior of the TA mode. In particular, when 
the nanotube is force-free the TA mode frequencies vary quadratically with wave number and have 
curvature (second derivative) proportional to the square-root of the nanotube's bending stiffness. When 
the nanotube has non-zero external force, the TA mode frequencies vary linearly with wave number 
and have slope proportional to the square-root of the axial force. In addition, we show that the out-of- 
plane shear deformation does not contribute to this mode and the unusual properties of the TA mode 
are exclusively due to bending. Therefore, the TA phonon curves - and associated transport properties 
- are not material properties but rather can be directly tuned by external loads. We illustrate this by 
showing the effect of the external axial load on the Density of States. We develop a simple mechanical 
model based on classical beam theory that captures the essential physics of the geometric nonlinearity 
in slender nanotubes that couples the axial load directly to the phonon curves. 

The paper is organized as follows. We first describe what is known in layered nanostructures in 
Section [2} We then present the closed- form calculations for atomic chains in Section |3j In Section]?} 
we present the results of numerical calculations for carbon nanotubes. Finally, in Section |5| we show 
that a simple model based on classical beam theory captures the key nonlinearity that governs the TA 
mode behavior. 

'The conclusions of this study appear mixed. In the authors' words: "We have shown that flexure modes, those with 
a; oc at long wave length, exist in carbon nanotubes. We also find that the flexure modes cross over to linear dispersion 
at very small values of wave vector. We are unsure whether this cross over is a feature of our choice of spring constants, or 
is an actual feature of nanotubes." 
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2 Layered crystals and atomic sheet structures 



Quadratic dependence on the wave vector for long wavelength TA phonons is well known in the context 
of layered crystals, such as graphite, and sheet structures, such as graphene. The first works in this 
area were by Komatsu [24ll25l (1951) and Lifshitz [i26il (1952). Their work showed that the TA phonon 
mode, with wave vector parallel to the layers, is a "bending mode" and has quadratic dispersion. Other 
works concerned with bending modes in layered crystals include [|27l|28l|29l[30l[3Tl. The dispersion 
relation results of these early investigations are summarized in the paper [|28l a^ 

Transverse in-plane mode: ul 

Longitudinal in-plane mode: U2 

Out-of-plane (bending) mode: a;| 

Here, 6*11,6*12,6*33 and 6*44 are the elastic constants of the crystal, p is the density, and ky are 
the components of the wave vector parallel to the layer planes and k^ is the component of the wave 
vector perpendicular to the layers. The quantity A characterizes the bending rigidity of layers and is 
determined mainly by intralayer forces. For layered crystals, the interlayer elastic constants 6*33 and 
C*44 are usually negligible compared to A. This illustrates that the predominant contribution to the 
bending mode dispersion relation is due to the bending rigidity A. 

It is important to note that Eqn. ([3]) holds for stress-free layers only. If the layer or sheet is subjected 
to non-zero in-plane stress then, according to Lifshitz, the dispersion of the bending mode is 

col = -p{kl + K) (4) 

where a is a measure of the stress in the layer. Thus, the TA bending mode dispersion relation of a 
stressed layer has a linear dependence on the long wavelength wave vector and has a slope proportional 
to the square -root of the stress. Making an analogy between layered crystals or sheets and nanotubes, 
one would then expect similar behavior of the TA phonon modes in nanotubes. In fact, this is exactly 
what we find in Section]?] from our numerical calculations. 



= '=^^^{kl + kl) + --fkl (1) 
= ^{kl + kl) + ^kl (2) 

= — {kl + kl) + —kl + ^'{kl + Kf- (3) 



3 Phonons in a Linear Atomic Chain 

In this section, we examine a linear atomic chain as a model system that permits closed-form calcu- 
lations. We show that the TA modes of a linear chain of atoms have linear dependence on the wave 
number when the chain is subjected to a non-zero axial force. Further, we show that the TA long 
wavelength phonon modes for di force-free chain correspond to bending, not shear. 

3.1 Transverse acoustic phonon modes 

For simplicity, we consider a linear chain of atoms which is free to deform in a two-dimensional space. 
In the configuration of interest the position vector of atom j is Xj = x"eQ, = ajei, where a is the 

^Note, these equations correspond to the resuhs of Lifshitz's work. Komatsu gives similar expressions except that he 
replaces the coefficients of the wave vector components by the squared sound velocities of graphite. 
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lattice spacing of the chain, ei is the axial direction of the chain, and a E {1,2}. A standard argument 
leads to the eigenvalue equation for the phonon frequencies 



(5) 



where m is the mass of each atom in the chain, k is the wave number (along the axial direction), cok 
is the phonon frequency, is the phonon polarization vector (eigenvector), and is the dynamical 
matrix given by 

Dfc = ^ Uoj exp [-ikx]] . (6) 



between 



Here, Upg is the Hessian of the chain's potential energy function = 0({xj}). That is, Upg = 
d'^(j)/dxpdxg. Due to the Euclidean invariance (objectivity) of the potential energy, we can — without 
loss of generality — take = 0({rpq}) to be a function of the pair-distances (r^y 
atoms in the chain [i32il . Then we find 

Ho,=EE 



dro, dxodxg 



EE 



i^O j^g 



droidrgj roiVg 
1 



(Xi - Xo) (g) (Xj - Xq 



(7) 



+ 6 ^ 



Xo) (g) (Xj - Xo 



5, 



droq rig 



(X„ - Xo 



(x„ - Xo 



In the standard basis, this becomes 



d'^ij) i j-q g 
droidrqj \i\ \j-q\ 

n A '9</' 1 n — A ^ JM_ J_ 

U "Og Z^i^o 9roi ^ "'^If droq \q\a. 



(8) 



From this form of the Hessian, it is clear that the polarization vector for the TA phonon modes of 
the linear atomic chain is vta = [0, 1]^- Thus for this case, the TA dispersion relation is found to be 



TA,fc 



~ y2 ^oj exp [-ikj 
j 

-E 



(9) 



oj \J\(^ 



-(1 — exp [—ikja]). 



Expanding the exponential in powers of (the small parameter) k and noting that dcjy/dvQj = d(l)/drQ_j, 
results in 



TA,k 



-E 



2m ^ droj 



Z-^ OArr 



24m dr, 



Oj 



k' + .... 



(10) 



Noting that \j\a = r^j and neglecting higher order terms we find 



The virial axial force in a chain of atoms at zero temperature can be shown 



to be F 



2L d^'^i^ ^ d^^^r Thus, the linear part of the dispersion relation (a; versus k) is 

proportional to the square-root of the axial force in the chain. Next, we show that the second term in 



( 1 1 ) is related to the bending stiffness of the chain. 



3.2 Energy change corresponding to shear deformation 

If 7 is the shear parameter in the infinite chain, the displacement of the atoms is m] = and m| = 7x]. 
The separation between the atoms in the deformed chain will be 

+ At,, = [rj + [u^ - u^f] = [rj + = r,, [l + 7^] - r,, ^ 

(12) 

Therefore the change in distance between the atoms is Avji ^ ^rj/ = ^ro(j_;). Expanding the 
potential energy density, (p = (j)/2Na, up to the first order in perturbation, gives an energy change of 
the form 



7' -7^ 



^ 2Na ^ dru ^ 



r 



7 d(p 



F72 



(13) 



where F is the virial axial force. As can be seen, the change in energy due to shear deformation is not 



related to the second term of ( 1 1 ). 



3.3 Energy change corresponding to bending deformation 

Now consider the chain subjected to bending. The atom with initial position = jaei will move to 
Yj = r[sm{aj /r)ei + cos(aj/r)e2], where r is the radius of curvature. The arc length does not change 
during the bending. If rji = (j — l)a is the separation between atoms j and / before the bending, after 
bending the separation will be 



r,i + Ar,,= [{y]-y]f + {y',-y]f] 



2r sin 



al 

1 — cos — cos 
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2] 1/2^ 
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This shows that 



2r 

24r2' 



4 



24^2- 



(14) 



(15) 



Expanding the potential energy density up to the first order in perturbation gives an energy change of 
the form 



M = Ar.il 
^ 2Na ^ dr. I ^ 



del) 



—y 



'(j) 3 —1 X - d(j) 3 



(16) 
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Hence, we find Z\0 ^ = Ak?, where 



A 



24a 



-1 



d(t) 

dror 



m 



r: 



,3 

On 



(17) 



is the bending stiffness and k. = 1/r is the bending curvature. Thus, we have shown that the second 
term on the right hand side of ( [TT] ) is related to the bending stiffness A. Consequently, the curvature of 
the force-free chain's TA dispersion relation is proportional to the square-root of its bending stiffness. 

4 Numerical Calculations for Carbon Nanotubes 

In this section, we use our recently developed symmetry-adapted phonon analysis method['23] to nu- 
merically demonstrate the connection between the TA dispersion relation and the axial force and bend- 
ing stiffness of a SWNT. All calculations are performed using the Tersoff interatomic potential! 35 1. 

Before we describe our results, we note that the paper [|36| shows, numerically, that bending modes 
with quadratic dispersion exist in quantum wires with rectangular cross section. The existence of only 
one continuous wave vector makes the wires similar to nanotubes. However, contrary to nanotubes, 
because of the rectangular cross section, the bending modes are no longer degenerate in wires. Our 
systematic symmetry-adapted phonon analysis method is able to deal with the degenerate behavior 
that occurs in nanotubes and clearly illustrates the quadratic dispersion of the TA phonon modes for a 
force-free tube. 

First, we investigate the effect of axial relaxation on the density of states (DOS) and the form of 
the dispersion relation for the bending modes. Then we numerically show that the slope and curvature 
of the TA modes correspond to the axial force and bending rigidity, respectively. 

Following the procedures outlined in our previous paper [|23l . the dispersion curves and DOS for 
a relaxed (force-free) (7,6) SWNT are plotted in Fig. [1} Adding axial force leads to one significant 
change of the dispersion curves and DOS. That is, the dispersion of low frequency TA modes changes 
from quadratic to linear in the long wavelength limit (as k — )■ 6), where 9 = O.SOTvr is the rotation 
angle of the (7, 6) nanotube screw generator. A consequence of the low frequency dispersion curve 
near A; = 6* is that the DOS is singular near cu = for the force-free tube, while conversely, it is 
almost constant for tubes with non-zero force. This effect is illustrated in Fig. |2} by plotting the long 
wavelength TA phonon dispersion relation near k = 9 for tubes subjected to four different values of 
axial stress. Here, we define the axial stress as the axial force divided by the force-free circumference 
of the tube. 

The figure shows clear evidence that the long wavelength TA phonon dispersion behavior is related 
to the axial stress in the SWNT. 

To show more conclusively that the slope and curvature of the long wavelength TA modes have 
the specific dependence postulated in this paper, we compute (1) the bending modulus of SWNTs as a 
function of tube diameter, (2) the slope of the long wavelength TA phonon mode dispersion relation as 
a function of axial stress and tube diameter, and (3) the curvature of the long wavelength TA phonon 
mode dispersion relations for a force-free SWNT as a function tube diameter. 

We compute the energy density versus curvature relation using zero-temperature objective molec- 
ular dynamics (based on [|22| which allows for the explicit control of curvature through a "bending 
group parameter") for armchair nanotubes with different diameters. In particular, we calculate the 
change in energy density, Z\0 = Z\0/(£P), where i is the objective cell length and P is the perimeter 
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Figure 1 : (Left) Dispersion curves and (right) DOS of zero-force axially relaxed (7,6) SWNT as de- 
scribed within the symmetry-adapted objective structures framework of paper [|23l 




Figure 2: The effect of axial stress on (left) the dispersion relation of the TA phonon modes and (right) 
the DOS of a (7,6) SWNT 
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Figure 3: The square root of bending stiffness (calculated by minimizing the energy of armchair nan- 
otubes) versus the diameter of the nanotube. The points are from calculations and the line is a best-fit 
to the points. 
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Figure 4: (a) Square of the slope of the TA mode vs. the axial stress, (b) Curvature of the TA mode vs. 
diameter for various nanotubes ranging from (6, 6) to (16, 16). The points are from calculations and 
the line is a best- fit to the points. 



of the nanotube. For each radius, we calculate the bending stiffness A according to 



Z\0 



1 



Ah" 



A 



2A(t) 



(18) 



Figure [s] shows the variation of versus the diameter. It is clear that the bending stiffness scales 
linearly with the square of the diameter. 

Next, Fig. ^a) shows the (signed) square of the TA modes' slope (near A; = ^^) as a function of 
axial stress for a variety of tube SWNT types. This slope was computed from our phonon dispersion 
calculations using a numerical differentiation algorithm. Finally, in Fig. ^b) we show the curvature, 
computed in a similar manner to the slope, for various nanotubes ranging from (6, 6) to (16, 16) as a 
function of the tube diameter. 

From these figures we can infer the following important results: 
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• Fig. shows that the slope of the TA mode curve does not depend on the specific nanotube. 
It is linearly related to the square root of the axial force. 

• For various nanotubes of different diamters and chiralities, Fig. |4]^b) shows that the curvature of 
the TA mode curve scales linearly with the tube diameter, and Fig. [3] in turn shows that square 
root of the bending stiffness scales linearly with the tube diameter. Together, these imply that 
the curvature of the TA mode curve is directly related to the square root of the bending stiffness. 

These plots therefore confirm our claim. 



5 Simple Continuum Model and Discussion 

First we consider the case of a rod that undergoes stretching deformations, but with a pre-stress. This 
gives us the expected conclusion that the pre-stress plays a role only in determining the equilibrium 
state around which linearization is conducted. Consider a rod of length L that is extend by an amount 
S and with a pre-stress (Tq. The energy of deformation is, to second order, given by 

6 1 / (5 \ ^ 
E^ao- + -ET(-] +... (19) 



L 2 

Here Et is the extensional modulus at the state of linearization. The natural frequencies from a normal 

1 

mode analysis will be proportional to the square root of the coefficient of the quadratic term, i.e. E^. 

Now, to understand the current setting, we examine a beam that undergoes transverse deflections. 
Consider a beam of length L, and impose a bending deformation that causes a relative deflection A 

between the ends of the beam. A simple accounting of the geometric nonlinearity shows that there 

1 

is an extensional strain of ^1 + (^)^ j ^ — 1 which is simply ^ (^)^ to leading order. The energy of 
deformation is, to second order, given by 

1 /Z\^ ^ 



E^-ao^-J +... (20) 

The natural frequencies from a normal mode analysis will be proportional to the square root of the 

1 

coefficient of the quadratic term, i.e. a§ . When this term is absent, higher-order effects such as 
bending control the curvature of the TA mode curve. 

Therefore, this gives the unusual direct dependence of the normal mode frequencies on the axial 
force. As the above simple analysis shows, the key to this feature is the geometric nonlinearity of a 
slender structure. 

As our calculations above have shown, the Density of States can also be tuned significantly by 
applying an axial force. Therefore, it is possible that properties such as heat transfer can be tuned by 
mechanical means. 
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